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Local skin friction coefficient 
Channel halfwidth 

Jacobian matrix of the generalized grid transformation 

Turbulent kinetic energy 

Length 

Maximum core Mach number in the Sajben diffuser 
Static pressure 
Exit static pressure 

Production rate of turbulent kinetic energy 
Inlet total pressure 

Pipe radius for axisymmetric pipe flows 

Pressure ratio for Sajben diffuser flows, Pe/Pt 

Reynolds number based on channel centerline velocity and halfwidth 

Reynolds number based on pipe centerline velocity and radius 

Reference Reynolds number, (p r u r l r )/[i r 

Reynolds number based on freestream velocity and distance along flat plate 
Reynolds number based on freestream velocity and momentum thickness 
Static temperature 
Inlet total temperature 
Physical time 
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Velocities in the Cartesian x and y directions for 2-D planar equations or in the 
cylindrical x and r directions for axisymmetric equations 
Velocities in the swirl direction in the axisymmetric equation 
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Cylindrical axial and radial coordinates 
Cartesian coordinates 

Minimum distance from the nearest solid wall 
Non-dimensional wall coordinates 
Boundary layer thickness 
Dissipation rate of turbulent kinetic energy 
Laminar and turbulent viscosity coefficients 
Static density 
Computational time 
Computational coordinates 
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Centerline flow properties 
Dimensional reference conditions 


INTRODUCTION 


Two-equation turbulence models are often used in conjunction with Reynolds-averaged Navier- 
Stokes computer codes to simulate turbulent flow fields, and there has been a considerable amount of 
research devoted to these models. However, previous turbulence research projects often concentrated 
on deriving and validating the turbulence models without devoting the considerable amount of efforts 
required to take a two-equation turbulence model from pure turbulence research to practical 
engineering calculations. As a result, there is a need to incorporate well proven two-equation 
turbulence models into a computer code that is user-oriented, well documented, and easy to use. It is 
hoped that this will permit fluid dynamicists to quickly and efficiently perform CFD calculations 
involving two-equation turbulence models for a wide variety of turbulent flow fields. In this paper, 
the implementation and validation of the Chien low Reynolds number k-e turbulence model (Chien, 
1982) in the two-dimensional/axisymmetric version of Proteus (Towne et al., 1990), a compressible 
Navier-Stokes computer code, will be described. 


Algebraic turbulence models are simple to use and require little computational time. As a result, a 
popular algebraic turbulence model, the Baldwin-Lomax turbulence model (Baldwin and Lomax, 
1978), is available in the Proteus code. However, it has generally been accepted that two-equation 
turbulence models are more general than algebraic models while still requiring a reasonable amount of 
computational time. In a review by Patel et al., 1985, the low Reynolds number two-equation models 
of Launder and Sharma, Chien, Lam and Bremhorst, and Wilcox and Rubesin are found to perform 
better than the others. For the current implementation, the Chien k-e model was chosen for the 
following reasons: 

• Its low-Reynolds number terms are approximated by simple algebraic relations. This reduces 
the computational time and improves the numerical stability of the k-e model in the near wall 
region. 

• Its boundaiy conditions for k and £ at the wall are simple. 

• It gives reasonable turbulent and mean flow results near a solid wall. 

• It has been used by other researchers in Navier-Stokes calculations with good results (Sahu 
and Danberg, 1986, Nichols, 1990). 

In this implementation, the Proteus philosophy of code readability and modularity has been 
strictly adhered to. To optimize performance on the Cray computers, the Fortran coding of the Chien 
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k-e model has been vectorized as much as possible. Additional logic has been provided to handle 
flows with multiple solid surfaces and axisymmetric flows with or without swirl. 


CHIEN k-e TURBULENCE EQUATIONS 


The two-dimensional/axisymmetric Chien k-e equations are transformed into generalized 
nonorthogonal body-fitted coordinates. For the 2-D case, the following generalized grid 
transformation is used to transform the k-e equations from physical (x,y,t) coordinates to 
computational coordinates: 

£ = £(x,y), V = Tl(x,y), and r = t (1) 


The transformed 2-D k-e equations can be written in vector notation as: 
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The turbulent viscosity is given by 

k 2 

P, = C M p- 
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The empirical constants are the same as those used in the original Chien formulations: 
G k = 1.0, cr E = 1.3 

C, = 1.35, C 2 = 1.8 


2 — 
l--e 36 


(14) 

(15) 

(16) 

(17) 

(18) 


C u = 0.09( 1 - e-* 0 " 5 ’* ), where R,=^~ 
v ’ pe 

Note that y n is the minimum distance to the nearest solid wall, and y + is computed from y n . The 
production of turbulent kinetic energy includes the full Boussinesq approximation for compressible 
flows. All of the above equations have been nondimensionalized using appropriate normalizing 
conditions (Anderson et al., 1984). The turbulent kinetic energy k has been nondimensionalized by 
u r 2 , and turbulent dissipation rate e by p r u r 4 /p r . 


Using the transformation of equation (1) with y replaced by r, the transformed axisymmetric k-e 
equations for the Chien model can be written using vector notation as 
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In the current implementation, the swirl velocity w is included in the axisymmetric formulation of 
the turbulent kinetic energy production term 


SOLUTION ALGORITHM FOR THE k-e EQUATIONS 

For code modularity in turbulence modeling, the set of k-e equations are lagged in time and solved 
separately from the Navier-Stokes equations using a finite difference technique. Like the Navier- 
Stokes equations in Proteus, the k-e equations are solved by marching in time using a coupled 
alternating direction implicit (ADI) solution procedure with the generalized first- or second-order time 
differencing scheme of Beam and Warming, 1978. The 2-D k-e equations are solved using the 
following sequence: 
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Sweep 2 (p direction) 
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An analogous solution procedure is used for the axisymmetric k-e equations. The parameters 0i 
and 02 determine the type of time differencing scheme used (Beam and Warming, 1978). To 
approximate spatial partial derivatives, a variable centered formula is used. Following Nichols, 1990, 

A A 

the spatial derivatives for the convective terms A, C, F., G c are approximated using first-order upwind 
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differencing. A first-order backward difference approximation is used for the terms with positive 
eigenvalues, and a first-order forward difference approximation is used for the terms with negative 
eigenvalues. Since the k-e solver uses the same time step as the Navier-Stokes solver, several time 
step options are available. In addition, the user can specify the number of k-e iterations per Navier- 
Stokes iteration and any k-e time step factor. 

The boundary condition types and boundary condition values for k and e can be specified at any 
boundary as a whole or at any point on any boundary. For easy modification and easy evaluation of 
complicated boundary condition expressions for k and e, the boundary conditions are treated 
explicitly. Spatially periodic boundary conditions for k and e are also available and implemented 
using a periodic block solver. Unlike the Proteus Navier-Stokes solver, artificial viscosities are not 
used in the k-e solver. Smoothing is provided by using first-order upwind differencing for the 
convective terms in the k-e equations as described above. No artificial bounds or constraints were 
placed on the computed k and e profiles. The only requirement imposed on the k and e profiles 
computed in this implementation is that they both have to be positive. 

For computations with the k-e turbulence model, initial values for k and e are required throughout 
the flow field to start the time marching procedure. The best choice for initial conditions for the k-e 
equations will vary from problem to problem, and the user must supply a subroutine that sets up the 
initial values for k and e for the time marching procedure. A default subroutine, subroutine KEINIT, 
is built into Proteus that computes the initial k and e values from an initial or restart mean flow field 
based on the assumption of local equilibrium (production equals dissipation), and variations of that 
scheme have been found to be useful in computing the k-e initial conditions for many turbulent flow 
calculations. 


In subroutine KEINIT, the production rate of turbulent kinetic energy P^ is first computed from 
the initial or restart mean flow field using equations (1 5)-( 17), assuming k is zero everywhere. Then 
the initial profile for the dissipation rate e is found using the following relation: 


e = 


Re r p 


(30) 


Next, equation (18) is used to compute the initial profile for the turbulent kinetic energy k, and the 
resulting k and e profiles are then used to start the time marching process for the k-e equations. 


validation 

To validate Proteus and the k-e model, laminar and turbulent computations have been done for 
several benchmark test cases: incompressible fully developed 2-D channel flow, fully developed 
axisymmetric pipe flow, boundary layer flow over a flat plate, and turbulent Sajben transonic diffuser 
flow. The laminar computations served as confidence checks for the Proteus code as well as test beds 
for the various combinations of mean flow initial and boundary conditions that would be useful for 
the turbulent cases. 

The turbulent computations were done with both the Baldwin-Lomax and the Chien k-e turbulence 
models. The same mean flow boundary conditions were used for laminar and turbulent flow 
calculations. Unless otherwise stated, the default artificial viscosity and numerical differencing 
schemes in Proteus were used, initial conditions for the k-e equations were generated using the 
subroutine KEINIT mentioned above, and convergence was assumed when the average absolute 
residual value for each of the mean flow equation became less than lxlO' 6 . If it was necessary to 
cluster grid points near a solid wall, the default gridpoint packing option in Proteus was used. This 
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option uses the Roberts transformation (Roberts, 1971) to pack points near specified boundaries. All 
of the computations were done on the Cray YMP at NASA-Lewis Research Center. 


FULLY DEVELOPED 2-D CHANNEL FLOW 

Figure la shows the grid geometry as well as the boundary conditions used for all of the fully 
developed 2-D channel flow test cases. The plates are separated by a distance 2H in the y direction 
and the mean flow is in the x direction. Because of symmetry, only half of the 2-D channel was 
modeled. 

L a min ar C h annel fiqw Te$ t Case 

Incompressible fully developed laminar flow in a 2-D channel with a centerline Mach number of 
0.05 was computed with no turbulence model. For this case, Ren is 100. A 5HxlH computational 
domain with 51x31 grid points was used. Uniform grid spacing was used in both the streamwise and 
normal directions. For mean flow initial conditions, the 2-D Poiseuille velocity profile was used. 
When a CFL number of 20 was used, the converged solution was obtained in 8710 iterations. As can 
be seen in figure 2a, the agreement between the Proteus result and the exact solution (Daily and 
Harleman, 1966) is excellent. For a given pressure drop, the centerline velocity computed by Proteus 
is 1.28% lower than the laminar exact solution, and this difference is observed to increase as the grid 
becomes coarser or the centerline Mach number is increased. This observation applies to all of the 
channel flow and pipe flow test cases discussed below. This computation required 2.20x10'^ 
sec/iteration/grid point. 

Turbulent Channel Flow Test Case 

The channel flow experiment conducted by Hussein and Reynolds, 1975 was selected for the 
turbulent fully developed 2-D channel flow test case. Identical flow conditions and grid point 
distribution were used for both the Baldwin-Lomax and the Chien k-e turbulence model. Ren used in 
this case is 32300. A 5HxlH computational domain with 31x61 grid points was used. Grid points 
were clustered in the normal direction near the channel wall with approximately 10 grid points located 
in the viscous sublayer, and the nearest grid point off the wall is at y + = 0.8. Uniform grid spacing 
was used in the streamwise direction. For mean flow initial conditions, the turbulent velocity profile 
derived by Musker, 1979 was used. When a CFL number of 40 was used, converged solutions were 
obtained in 43380 iterations for the Baldwin-Lomax model and 44860 iterations for the Chien k-e 
model. The large number of iterations required for convergence is caused by using a compressible 
code to compute an incompressible flow. 

Results in figures 2b, 2c, and 2e show that the Chien k-e model gave better predictions for the 
internal 2-D channel flow than the Baldwin-Lomax model. In this test case, the centerline velocities 
computed by Proteus using the Baldwin-Lomax model and the Chien k-e model are 0.34% and 3.3% 
lower than the experimental value, respectively. From figure 2c, it can be seen that the higher 
centerline velocity computed by the Baldwin-Lomax model is caused by the prominent wake-like 
behavior of the u + profile near the centerline. 

In figure 2d, the turbulent kinetic energy profile computed by the Chien k-e model is compared 
with the experimental data from Laufer, 1949 at Ren = 30800. The computed outer turbulent kinetic 
energy profile is slightly higher than the experimental data, but the trend is well predicted. While the 
near wall data published by Laufer was incomplete, the Chien k-e model predicted a prominent peak 
of turbulent kinetic energy in the near wall region as expected. Figure 2e shows that while the eddy 
viscosity profile predicted by the Chien k-e model gives better agreement with the experimental data 
than the Baldwin-Lomax model, it is significantly larger than the experimental eddy viscosity profile 
near the channel center. This discrepancy is also observed by Myong and Kasagi, 1990 in their 
evaluation of the Launder-Sharma and the Lam-Bremhorst low Reynolds number k-e models. The 
turbulent computations with the Baldwin-Lomax model and the Chien k-e model required 2.36x1 0' 5 
sec/iteration/grid point and 3.21xl0' 5 sec/iteration/grid point, respectively. 
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FULLY DEVELOPED AXISYMMETRIC PIPE FLOW 


Figure lb shows the grid geometry as well as the boundary conditions used for all of the fully 
developed axisymmetric pipe flow test cases. The pipe radius is R, and the mean flow is in the x 
direction. Because of symmetry, only half of the pipe was modeled. Both the Navier-Stokes and the 
k-e solvers were used in axisymmetric mode. 

Laminar Pipe Flow Test Case 

Incompressible fully developed laminar flow in an axisymmetric pipe with a centerline Mach 
number of 0.05 was computed with no turbulence model. For this case, ReR is 100. A 5RxlR 
computational domain with 51x31 grid points was used. Uniform grid spacing was used in both the 
axial and radial directions. For mean flow initial conditions, the Poiseuille velocity profile was used. 
When a CFL number of 20 was used, the converged solution was obtained in 3500 iterations. As can 
be seen in figure 3a, the agreement between the Proteus result and the exact solution (Daily and 
Harleman, 1966) is excellent for both the velocity profile and pressure drop. For a given pressure 
drop, the centerline velocity computed by Proteus is 0.58% lower than the laminar exact solution. 

This computation required 2.49xl0" 5 sec/iteration/grid point. 

Turbulent Pipe Flow Test Case 

The pipe flow experiment conducted by Laufer, 1952 was selected for the fully developed 
axisymmetric pipe flow test case. Identical flow conditions and grid point distribution were used with 
both the Baldwin-Lomax and the Chien k-e turbulence model. The ReR used in this case is 250000. 

A 5RxlR computational domain with 21x61 was used. Grid points were clustered in the radial 
direction near the pipe wall with approximately 8 grid points located in the viscous sublayer, and the 
nearest grid point off the wall is at y+ = 1.2. Uniform grid spacing was used in the axial direction. 
For k-e initial conditions, the default subroutine KEINIT described above was used. When a CFL 
number of 40 was used, converged solutions was obtained in 17280 iterations for the Baldwin- 
Lomax model and 13180 iterations for the Chien k-e model. 

The comparisons in figures 3b-3e show good agreements between experimental and 
computational results for both the Baldwin-Lomax model and the Chien k-e model. As in the 
turbulent channel flow test case, using the Chien k-e model resulted in better mean velocity profiles 
than the Baldwin-Lomax model (figures 3b-3c). In figure 3d, the turbulent kinetic energy profile 
computed by the Chien k-e model is slightly higher compared with the experimental data from Laufer, 
1952, but the peak near the wall and the general trend are well predicted. Figure 3e shows that the 
Reynolds stress profiles predicted by the Chien k-e and the Baldwin-Lomax models agree exactly 
with the experimental data. In this figure, the Baldwin-Lomax Reynolds stress profile is exactly on 
top of the Chien k-e profile. 

In this test case, the centerline velocities computed by Proteus using the Baldwin-Lomax model 
and the Chien k-e model are 2.09% higher and 2.34% lower than the experimental value, 
respectively. As in the turbulent channel flow test case, the Baldwin-Lomax model predicted a 
somewhat higher centerline velocity than the Chien k-e model. The turbulent computations with the 
Baldwin-Lomax model and the Chien k-e model required 3.71xl0' 5 sec/iteration/grid point and 
4.36x1 0' 5 sec/iteration/grid point, respectively. 


2-D BOUNDARY LAYER FLOW 

Figure lc shows the grid geometry as well as the boundary conditions used for all of the 2-D 
boundary layer test cases. Incompressible laminar and turbulent flows over a flat plate at zero 
pressure gradient with a freestream Mach number of 0.1 was computed with and without the 
turbulence models. The mean flow is in the x direction. For the upstream boundary, the Blasius 


8 


velocity profiles were used to specify the u and v boundary conditions for both the laminar and 
turbulent flow cases. To keep approximately the same number of grid points inside the boundary 
layer at any given x station, the height of the grid system is increased with the distance downstream of 
the plate. The grid points are clustered near the solid wall so that, at any given streamwise location, 
approximately 90 percent of the grid points are located inside the boundary layer. Because of this 
special grid treatment, the default Proteus grid packing scheme was not used. Instead, an algebraic 
grid packing scheme described in Hoffmann, 1989 was employed. 

Laminar Boundary Laver Flow Test Case 

Incompressible laminar boundary layer flow over a flat plate was computed. The computational 
domain ranges from Re x = 1000 at the upstream boundary to Re x = 100000 at the downstream 
boundary with 51 grid points in the x direction and 31 grid points in the normal direction. Grid 
points were clustered in the normal direction near the plate surface. Uniform grid spacing was used 
in the streamwise direction. For mean flow initial conditions, the Blasius velocity profiles were used. 
When a CFL number of 20 was used, convergence was reached in 1860 iterations. As can be seen in 
figures 4a and 4b, the agreement between the Proteus result and the exact solution (Daily and 
Harleman, 1966) is excellent for both the velocity and local skin friction profiles. This computation 
required 2.07 lxlO' 5 sec/iteration/grid point. 

Turbulent Boundary Laver Flow Test Case 

The incompressible turbulent boundary layer experiment conducted by Klebanoff, 1953 was 
selected for the turbulent boundary layer flow test case. Identical flow conditions and grid point 
distribution were used with both the Baldwin-Lomax and the Chien k-e turbulence models. The 
computational domain ranges from Re x = l.OxlO 4 at the upstream boundary to Re x = lO.OxlO 6 at the 
downstream boundary with 81 grid points in the x direction and 51 grid points in the normal 
direction. Grid points were clustered in the normal direction near the plate surface such that, at Ree = 
7700, approximately 10 grid points were located in the viscous sublayer with the nearest grid point 
off the wall at y + = 0.5. Uniform grid spacing was used in the streamwise direction. In the 
calculation using the Baldwin-Lomax model, the Blasius velocity profiles were used for the mean 
flow initial conditions. The CFL values were 5 for the first 300 CPU seconds and 20 thereafter. The 
converged solution was obtained in 12010 iterations. The computation using the Chien k-e model can 
also be started from the Blasius velocity profiles, but to save Cray CPU time, the Chien k-e 
calculation was started from the converged Baldwin-Lomax solution. The CFL values used in the 
Chien k-e computation were 2 for the first 300 CPU seconds and 10 thereafter. An additional 1 1920 
iterations were required before all of the residual values for each of the mean flow equations leveled 
out. It should be noted that the number of iterations required for convergence can be significantly 
reduced if higher freestream Mach numbers were used. For instance, the number of iterations 
required for convergence was cut in half when the freestream Mach number was increased to 0.3 with 
virtually no change in the computed results. 

Results in figures 4c-4d show that, although there are minor differences between the Baldwin- 
Lomax and the Chien k-e computed velocity profiles, both of those models predicted velocity profiles 
that are in excellent agreement with the experimental data. 

In figure 4e, the turbulent kinetic energy profile computed by the Chien k-e model is compared 
with the experimental data from Klebanoff, 1953. The agreement is very good. In figure 4f, the 
Chien k-e model is seen to give better prediction of the local skin friction coefficient profile than the 
Baldwin-Lomax model. The flat plate computations with the Baldwin-Lomax model and the Chien k- 
e model required 1.94xl0' 5 sec/iteration/grid point and 2.73xl0 -5 sec/iteration/grid point, 
respectively. 
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SAJBEN TRANSONIC DIFFUSER FLOW 


Figure Id shows the grid geometry as well as the boundary conditions used for all of the Sajben 
transonic diffuser test cases. Turbulent flow was computed in the converging-diverging duct. The 
mean flow is in the x direction. Given a fixed inlet total pressure, there can be a strong shock, weak 
shock, or no shock at all in the diffusing section after the throat, depending on the pressure at the 
diffuser exit. The weak shock case has been used as a test case in the 2-D Proteus code User's Guide 
(Towne et al., 1990). In this report, all three cases were computed using both the Baldwin-Lomax 
model with the Launder-Priddin mixing length formula in the inner region and the Chien k-e model. 

A computational domain of 81 grid points in the streamwise direction and 51 grid points in the normal 
direction was used. To capture the shock in transonic flow cases, more grid points were positioned 
in the streamwise direction in the diffusing section after the throat. Grid points were also packed in 
the normal direction near the two solid walls to resolve the viscous regions with the nearest grid point 
off the wall at y + = 1.5 - 4.0. The same grid distribution was used in all three cases. Since second- 
order time differencing was use in all of the calculations, no residuals was computed. Convergence 
was assumed when the pressure distributions remained essentially constant for a given fixed pressure 
drop. 

Initially, the Baldwin-Lomax model was used to start the calculation for all three cases. The initial 
conditions were zero velocity and constant pressure and temperature everywhere in the flow field. 

For the transonic flow cases, the exit pressure was gradually lowered to R = 0. 1338 to establish 
supersonic flow in the diffusing section after the throat. Then the exit pressure was gradually raised 
to the appropriate pressure ratio, R (R = 0.72 for strong shock and R = 0.82 for weak shock) to 
establish the shock and iterated there using either the Baldwin-Lomax or the Chien k-e turbulence 
models. For the subsonic flow case with no shocks, starting from zero initial conditions everywhere, 
the exit pressure was gradually lowered to R = 0.862 using the Baldwin-Lomax model and then 
iterated there using either of the turbulence models. 

Figures 5a-5f compare the pressure distributions computed by Proteus using the Baldwin-Lomax 
and the Chien k-e models with the experimental data used by Hsieh et al., 1987. For the no shock 
case (R = 0.862), the Chien k-e results agree well with the experimental data while the Baldwin- 
Lomax model gave pressures that are too low at the throat. Both turbulence models give good results 
for the weak shock case (R = 0.82). For the strong shock case (R = 0.72), the Chien k-e model does 
a better job of predicting the pressure profile right after the shock. However, both turbulence models 
predicted pressures that are too high after the shock. This is to be expected, since experimental data 
published by Hsieh et al., 1987, Bogar et al., 1983 and Salmon et al., 1983 showed that flow 
separation does occur at the upper wall for the strong shock case, and no special modifications were 
made to either of the models to compute separated flow. Note that the shock locations for the 
transonic cases could be fitted better if the grid points were clustered at the shock locations. But since 
this was a shock capturing analysis, identical grids were used in all cases. 

Figures 6a-6f plot the computed Mach number contours for the no shock, weak shock, and strong 
shock cases, respectively. The Mach contours for the no shock and the weak shock cases computed 
using the Baldwin-Lomax model look similar to those computed using the Chien k-e model, but the 
Baldwin-Lomax model predicted a maximum core Mach number that is too high for both cases. For 
the strong shock case, the Baldwin-Lomax model incorrecdy predicted bottom wall flow separation 
after the shock, while the Chien k-e calculation and experimental results showed that separation 
occured on the top wall. Using either of the turbulence models, Proteus computed the critical mass 
flow rate for the transonic diffuser flow cases to be about 1-3% lower than the one-dimensional 
isentropic equation. For the Sajben diffuser, calculations using the Baldwin-Lomax model and the 
Chien k-e model required 4.00x1 0' 5 sec/iteration/grid point and 4.24x1 0' 5 sec/iteration/grid point, 
respectively. 
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SUMMARY OF VALIDATION CASES 


From the above validation cases, one can see that the Chien k-e model consistently predicted the 
correct flow physics for all of the flows under consideration, whereas the Baldwin-Lomax model was 
having problems with the internal flow cases. The Chien k-e model required about 41% more CPU 
time than the Balwin-Lomax model for simpler flows such as the flat plate and only about 6% more 
CPU time for more complicated flows such as the Sajben diffuser. For more complicated flows, 
especially flows with multiple solid walls, the difference in CPU times required by the Chien k-e 
model and the Baldwin-Lomax model is small. This is because the Baldwin-Lomax model is 
burdened with repetitious calculations, additional logic, and complicated averaging processes in order 
to come up with a reasonable turbulent viscosity profile for flows with complex geometries. Whereas 
in the Chien k-e model, more complex flow cases use exactly the same amount of calculations as the 
simpler flow cases with only changes in boundary conditions. 


CONCLUSIONS 

The Chien k-e turbulence model has been successfully implemented into the 2-D version of 
Proteus, a general purpose Navier-Stokes code. Validation test results showed that the Chien k-e 
model gives good results for a wider variety of flows than the Baldwin-Lomax model in Proteus. 
Also, the above test cases showed that the CPU time required for Proteus calculations using the Chien 
k-e model is only slightly more than the Baldwin-Lomax model for real life flow cases with 
complicated geometries. 
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FIGURES 
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Fig. la Grid geometry and boundary conditions for fully developed 2-D channel flow 
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Fig. 1b Grid geometry and boundary conditions for fully developed axisymmetric pipe flow 
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Fig. 1c Grid geometry and boundary conditions for boundary layer flow over a flat plate 
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Fig. Id Grid geometry and boundary conditions for Sajben diffuser flow 
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Fig. 2a Fully developed laminar channel flow 
velocity profile 



Fig. 2c Fully developed turbulent channel flow u + profile 
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Fig. 2b Fully developed turbulent channel flow 
velocity profile 
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Fig. 2e Fully developed turbulent channel flow eddy 
viscosity profile 
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Fig. 3a Fully developed laminar pipe flow velocity profile 



Fig. 3c Fully developed turbulent pipe flow u + profile 
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Fig. 3e Fully developed turbulent pipe flow Reynolds 
stress profile 



Fig. 3b Fully developed turbulent pipe flow velocity profile 



Fig. 3d Fully developed turbulent pipe flow 
kinetic energy profile 


16 



Fig. 4a Laminar boundary layer velocity profile 



Fig. 4b Laminar boundary layer local skin friction 
coefficient distribution 



U/L L 

Fig. 4c Turbulent boundary layer velocity profile 
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Fig. 4d Turbulent boundary layer u + profile 



Fig. 4e Turbulent boundary layer kinetic energy profile 


Fig. 4f Turbulent boundary layer local skin friction 
coefficient distribution 
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p/p p/p p/p 



Fig. 5a Sajben diffuser top wall pressure distribution, 
no shock case 



Fig. 5b Sajben diffuser lower wall pressure distribution, 
no shock case 




Fig. 5c Sajben diffuser top wall pressure distribution, Fig. 5d Sajben diffuser lower wall pressure distribution, 

weak shock case weak shock case 



x/H 

Fig. 5e Sajben diffuser top wall pressure distribution, 
strong shock case 



Fig. 5f Sajben diffuser lower wall pressure distribution, 
strong shock case 


18 



19 


Fig. 6 Computed Mach number contours for the Sajben diffuser 
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